Orbital control of noncoUinear magnetic order in nickelate heterostructures 
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We have used resonant x-ray diffraction to develop a detailed description of antiferromagnetic ordering in 
epitaxial superlattices based on two-unit-cell thick layers of the strongly correlated metal LaNiOg. We also 
report reference experiments on thin films of PrNiOg and NdNiOg. The resulting data indicate a spiral state 
whose polarization plane can be controlled by adjusting the Ni d -orbital occupation via two independent 
mechanisms: epitaxial strain and quantum confinement of the valence electrons. The data are discussed in 
the light of recent theoretical predictions. 
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The theoretical prediction and experimental control of col- 
lective electronic ordering phenomena in metals are among 
the greatest challenges of current solid-state physics. Re- 
cent advances in the synthesis of epitaxial metal-oxide het- 
erostructures have provided new^ opportunities for the ex- 
ploration and manipulation of the phase behavior of metal- 
lic conductors [1, 2]. How^ever, detecting the subtle spin, 
charge, and orbitally modulated structures resulting from 
Fermi surface instabilities in the atomically thin layers of 
such devices remains challenging. Most experiments to-date 
have focused on spatially uniform order parameters (such 
as the uniform magnetization) that modify the macroscopic 
properties, and on corresponding layer-to-layer variations 
w^hich can be determined by neutron or x-ray reflectome- 
try. In contrast, little experimental information is currently 
available about ordering patterns within the heterostructure 
layers, despite many predictions in this regard. 

We have used resonant x-ray diffraction (RXD) to de- 
velop a detailed, microscopic description of the magnetic 
ordering patterns in epitaxial superlattices (SLs) of metal- 
He perovskites LaNiOg (LNO)and thin films of RNiOg (R = 
Nd and Pr) . Bulk analogs of these materials either remain 
paramagnetic at all temperatures (R = La), or become Mott- 
insulating and antiferromagnetic at low temperatures (R / 
La) [3, 4]. Recent experimental [5-10] and theoretical [11- 
19] work has focused on the phase behavior of RNi03 -based 
thin films and heterostructures. In particular, muon spin ro- 
tation (juSR) experiments yielded evidence of antiferromag- 
netism in SLs with LNO layers of thickness below three unit 
cells (u.c.) sandwiched between insulating blocking layers, 
whereas SLs with thicker LNO layers remain paramagnetic 
[5]. Exchange-bias effects revealed by magnetometric mea- 
surements on heterostructures of LNO and ferromagnetic 
LaMn03 also indicated antiferromagnetic or spin-glass order 
in LNO [6]. These results have established the nickelates as 
a promising platform for phase control of metals in proxim- 
ity to a Mott transition, and they have stimulated detailed 
theoretical predictions for the magnetic ordering patterns 
in nickelate heterostructures and SLs [11, 12, 18, 19]. In 
analogy to prior work on bulk-like NdNi03 (NNO) films 



[20], we have used RXD to show that the Ni moments in 
the SLs form a spiral state. We further demonstrate control 
of the polarization plane of the spiral by manipulating the 
Ni d -orbital occupation via epitaxial strain and quantum 
confinement. Our RXD data enable detailed tests of theo- 
retical work on the nickelates, and provide the foundation 
for the integration of "orbitally engineered" nickelates into 
spintronic devices. 

Superlattices of LNO and insulating, nonmagnetic LaAl03 
(LAO) or DySc03 (DSO) [21] were grown using pulsed laser 
deposition on various [001] -oriented substrates: La(Sr)Al04 
(LSAO), [LaAl03]o.3[Sr2AlTa06]o.7 (LSAT), SrTi03 (STO), 
and DSO with in-plane lattice constants 3.756, 3.87, 3.905, 
and 3.95 A, respectively. The substrates and blocking layers 
allowed us to impose variable levels of compressive or ten- 
sile strain on LNO, whose bulk lattice constant is a = 3.84 
A. Data on NNO and PrNi03 (PNO) films with a = 3.79 and 
3.84 A, respectively, were taken for comparison. We discuss 




Fig. 1 (color online). Left panel: Schematic diagram of the mag- 
netic structure of a AT = 2 LNO-RXO SL grown on tensile strain. 
The Ni moments are illustrated by red arrows. Also shown are the 
electron density distributions in the Ni dx2_y2 (opaque) and d^z^-r^ 
(semitransparent) orbitals. Right panel: samples were mounted 
on a tilted wedge to access the [111] direction of the pseudocubic 
perovskite structure in the horizontal scattering geometry. Also 
shown are the wave vectors of the incoming and outgoing photons 
(green arrows), the corresponding momentum transfer q, the in- 
coming photon polarization s (black arrows), and the azimuthal 
angle ip. 
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Fig. 2 (color online), (a) Scans around qgow = (t? Tj ^) for SLs with AT = 2 consecutive LNO unit cells show a magnetic Bragg reflection, 
while those with AT = 3, 4 do not. The inset shows a reciprocal-space map of the scattered intensity from the N=2 SL. (b) Comparison 
between the x-ray absorption spectrum (XAS) and the photon-energy dependence of the magnetic Bragg intensity at qsow close to the Ni 
L3 (2p -^ 3d) edge, which shows almost identical line shapes for two azimuthal angles ijj = 0,30°. (c) Temperature dependence of the 
magnetic Bragg intensity at qg^w ^^ LNO -based SLs with AT = 2 (red symbols), compared to representative dc electrical conductance 
measurements (taken by a standard four-probe method) for AT = 2 and AT = 4 SLs. 



structural data in terms of the pseudocubic unit cell of the 
perovskite lattice. The momentum transfer q = (H,K,L) 
is also indexed in these units. As demonstrated previously 
[5, 7], the layer sequence of our SLs is defined w^ith u.c. 
precision. We consider SLs w^ith N u.c. LNO - N u.c. RXO 
(A/ = 2,3,4) structures, labeled N//N) in the foUow^ing. 

The RXD experiments w^ere performed at the BESSY-II 
undulator beam line UE46-PGM1 using variable linearly po- 
larized photons in a horizontal scattering geometry (Fig. 1). 
A three-circle UHV diffractometer w^as equipped w^ith a 
continuous -flow^ He cryostat. To access the i^AA) reflec- 
tion, the [001] -oriented samples v\^ere mounted on copper 
wedges with a 55° tilt. The azimuth value 1/^ = 0° is assigned 
when the crystal vectors [1, 1, 1] and [1, 1, 0] span the scat- 
tering plane, and positive rotation around q is left-handed. 
This geometry precluded full 1/^ -scans, as there are two angu- 
lar ranges where the incident or outgoing beams are below 
the sample horizon. In our definition of ip, these occur 
around ip = 270° and 90°, respectively. To estimate the val- 
ues of the polarization-dependent scattering cross section at 
these azimuth positions, we used a non-tilted (asymmetric) 
scattering geometry to access part of the (^, ^, L) scattering 
rod. 

Figure 2(a) shows selected reciprocal-space scans on 
(2//2), (3//3), and (4//4) SLs taken under resonant con- 
ditions, that is, with the photon energy tuned to the Ni 
L3 edge. (2//2) SLs of all compositions grown on all sub- 
strates investigated here show resonant Bragg reflections 
at q = (^, ^,L), whereas samples with 3 or 4 consecutive 
u.c. of LNO do not. Comparison of the energy dependence 
of the diffracted intensity at this q with x-ray absorption 
spectroscopy (XAS) data (Fig. 2(b)) confirms the resonance 
at the Ni L3 edge, and the polarization dependence of the 
scattering cross section (see below) confirms the magnetic 
origin of the Bragg reflections. The RXD data thus yield 
the ordering vector of the antiferromagnetic state identified 



by juSR [5]. The observation of magnetic Bragg reflections 
implies staggered ordering of the Ni moments in the SL 
plane analogous to the one characterizing antiferromag- 
netic order in bulk nickelate perovskites [4, 20] and rules 
out a spin glass state that has been considered in related 
SLs [6]. The observation of a scattering rod perpendicu- 
lar to the SL plane (inset in Fig. 2(a)) is consistent with 
quasi-two-dimensional order, but the accessible range of the 
corresponding momentum transfer, L, is not sufficient to 
determine the spin-spin correlation length in this direction. 
Figure 2 (c) demonstrates that the onset of the magnetic 
Bragg intensity as a function of temperature coincides with 
the maximum observed in dc conductance measurements of 
N = 2 SLs. This observation agrees with recent predictions 
according to which the magnetic order in LNO -based SLs 
is a consequence of spin density wave (SDW) formation 
[11, 12]. We thus refer to the magnetic ordering vector as 
^SDw ' Further work is required to elucidate whether this 
transition is accompanied by a modification of the charge 
density [22]. 

In order to determine the SDW polarization, we systemat- 
ically varied the incoming and outgoing photon polarization 
vectors, s and s\ relative to the electronic magnetic mo- 
ments by scanning the azimuthal angle xjj around qsDw 
(Fig. 1). The intensity of a magnetic Bragg reflection is 
given by [23] 



hi=\J]e^''-^^e;'Fj(E)-e, 



(1) 



with k,l ^ {a,n}. Figure 2(b) shows that the energy profiles 
of the Bragg reflections measured at different azimuthal 
angles ip are identical. This confirms that the Bragg intensity 
at qsDw arises from a single set of scattering tensors Fj(E) 
at lattice site j, which (together with the corresponding 
phase factors) encode the magnetic ordering pattern. We 
measured the scattered intensity for a- and 7i-polarized 
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Fig. 3 (color online), (a)-(c) Dependence of the polarization dependent scattering intensity k/ct on the azimuthal angle ip for LNO-based 
SLs and a PNO film on LSAT (tensile strain), (b) a NNO film on LSAO (compressive strain, unconfined), and (c), a LNO-LAO SL on 
LSAO (compressive strain, confined). The solid and dashed lines represent the results of fits to non-collinear (NCM) and collinear (CM) 
magnetic structure models described in the text, respectively. In the shaded areas around ijj = 90°, measurements were not possible 
in the tilted scattering geometry of Fig. 1, and an alternative geometry was used (see text), (d)-(f) Sketches of the spin directions in 
the two-sublattice magnetic structure resulting from the fits. The color-coding refers to panel (a)-(c) [27]. (g)-(i) Linear dichroism 
obtained from XAS measurements around the Ni L2 edge with incoming light polarization in and out of the substrate plane. The ratio X 
of orbital occupations derived from the sum rule (see text) is given in the legend and schematically depicted by color-coded solid and 
semitransparent plots of the corresponding electron density distribution. 



incident photons (Fig. 1), so that the tw^o channels n = 
I^^ + 1^^ and a = I^j^ were distinguished. We present the 
measured data as the ratio n/a, which is not influenced by 
the shape of the sample and the orientation of its surfaces 
relative to the incoming and outgoing x-ray beams, and can 
thus be directly compared to model calculations. 

Figure 3 displays the polarization-dependent data for 
different LNO-based SLs with N = 2 and films exhibiting 
SDW order, grown on substrates that impose either tensile 
(Fig. 3(a)) or compressive (Figs. 3(b,c)) strain. The blue 
line in Fig. 3 (c) represents the result of a calculation for 
the magnetic structure of bulk RNiO^ with R^La [4, 20], 
which comprises two sets of spirals polarized perpendicular 
to the propagation vector qsnw (blue arrows in Figs. 3(f)). 
Whereas the magnetic structure of a LNO-LAO SL under 
compressive strain is compatible with the one in the bulk 
(Fig. 3(c)), the data sets on all samples grown under tensile- 
strain conditions (Fig. 3(a)), and the one on a NNO film 
grown under compressive strain (Fig. 3(b)), indicate dis- 



tinctly different magnetic structures. 

In analyzing the azimuthal scans of Fig. 3, we consid- 
ered both collinear antiferromagnetic structures with spatial 
variation of the moment amplitude and spiral structures 
analogous to the one in bulk RNi03 with identical ampli- 
tude on every lattice site (see Fig. 3(a,b)). Model calcula- 
tions [11, 12] show that the relative stability of these two 
structures depends on factors that are difficult to compute 
from first principles, such as the on-site correlation strength. 
The comprehensive data sets displayed in Fig. 3 allow an 
experimental test of these predictions. Whereas collinear 
structures turned out to be incompatible with the data (see 
dashed lines in Fig. 3(a,b)), calculations based on spiral 
SDWs yield excellent descriptions of all three distinct data 
sets (solid lines in Figs. 3(a,b))). These structures can be 
derived from the structure of bulk RNi03 by adjusting the 
direction of the moments in the two sublattices. In the best 
fits for all samples under tensile strain, the moments are 
symmetrically tilted from the [001] -axis by 28 =b 2° (red 



arrows in Fig. 3(d)), and remain coplanar with those of the 
bulk structure [27]. The magnetic structure of the NNO 
sample under compressive strain, on the other hand, com- 
prises moments along [110] and [110] (green arrows in 
Fig. 3(e)). Note the striking difference between the SDW 
polarization of this sample and the one of the LNO-LAO SL 
shown in Fig. 3(f), which was grown on the same substrate 
(LSAO) and exhibits the same in-plane lattice parameter. 

In order to uncover the origin of the surprising variability 
of the SDW polarization, we used x-ray absorption spec- 
troscopy with linearly polarized x-rays near the Ni L2 edge 
to determine the relative occupation of the Ni d-orbitals, 
which controls the magneto-crystalline anisotropy via the 
spin-orbit coupling [24, 25]. Figs. 3(g-i) display the differ- 
ence of the absorption spectra for photons polarized parallel 
and perpendicular to the substrate surface for three repre- 
sentative samples. Taking advantage of the sum rule for lin- 
ear dichroism, we have converted the energy-integrals into 
the ratio of e^ hole occupation numbers X = n2z^_r2/n^2_ 2 
[26]. The SLs and films grown under tensile strain show 
1.03 < X < 1.14, corresponding to an enhanced electron 
occupation of the (1^2 _y2 orbital (Fig. 3(g)). Since orbital 
moments in this situation will point along z, the experimen- 
tally observed canting of the spin moments towards this 
direction (Fig. 3(d)) is a natural consequence of the intra- 
atomic spin-orbit coupling. Conversely, the preferential 
^3z2-r2 occupation (X = 0.97) found in the compressively 
strained NNO film (Fig. 3(h)) accounts for the observed spin 
polarization in the x,y-plane (Fig. 3(e)). In the compres- 
sively strained LNO-LAO SL, our data show that the equal 
population of d^2_y2 and ^322. ^2 found in bulk nickelates is 
restored (Fig. 3(i)), consistent with prior reports on related 
SLs [9] and with the bulk-like magnetic structure of this sam- 
ple (Fig. 3(f)). This requires a mechanism that counteracts 
the effect of compressive strain observed in the NNO film 
on the same substrate. Such a mechanism has been identi- 
fied in model calculations for LNO-based SLs, [13, 14, 16] 
which indicate that the confinement of the LNO conduction 
electrons to the x,y -plane by the insulating blocking layers 
of the SL stabilizes the d^2_y2 orbital relative to 0^322 _ ^2. In 
contrast to the long-range effect of epitaxial strain, recent 
results from resonant x-ray reflectometry [7, 28] indicate 
that the orbital polarization is spatially modulated due to 
the dimensional confinement which predominantly acts on 
the interface layers, thus affording an independent means 
of controlling the SDW polarization. 

In summary, our RXD data have allowed us to develop a 
detailed picture of a complex in-plane magnetic superstruc- 
ture in an oxide superlattice. Along with recent neutron 
diffraction work on manganate superlattices [29], this rep- 
resents one of few cases where such a comprehensive de- 
scription has proven possible. Note that neutron diffraction 
is not applicable to systems with few magnetically ordered 
atomic layers, such as our SLs and most other oxide het- 
erostructures. Our demonstration that RXD experiments 
are feasible thus indicates new perspectives for this field of 
research. We have further demonstrated understanding and 



control of the noncollinear spin polarization by virtue of two 
independent mechanisms, epitaxial strain and quantum con- 
finement. Since the nickelates can be epitaxially integrated 
with other metallic oxides exhibiting ferromagnetic or su- 
perconducting order, the control options we have identified 
provide interesting new perspectives for spintronic devices. 
In view of the substantial conductance in the SDW state, 
these may eventually include devices with electronically 
active antiferromagnetic elements [30] . 
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